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Interacting two-component Fermi gases loaded in a one-dimensional {ID) lattice and subject to 
harmonic trapping exhibit intriguing compound phases in which fluid regions coexist with local Mott- 
insulator and/or band-insulator regions. Motivated by experiments on cold atoms inside disordered 
optical lattices, we present a theoretical study of the efltects of a random potential on these ground- 
state phases. Within a density-functional scheme we show that disorder has two main effects: (i) it 
destroys the local insulating regions if it is sufficiently strong compared with the on-site atom-atom 
repulsion, and (ii) it induces an anomaly in the compressibility at low density from quenching of 
percolation. 
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Introduction — The interplay between interactions and 
disorder in quantum many-body systems is an area of 
long-standing interest. For instance, both long-ranged 
Coulomb interactions and disorder from various mecha- 
nisms are believed to play an important role in the metal- 
insulator transition (MIT) in the two-dimensional (21?) 
electron liquid in semiconductor heterostructures"'^. Dis- 
order and interactions affect not only transport proper- 
ties of the 2D liquid, but also thermodynamic quantities 
such as the compressibility^'^ and the spin susceptibil- 
itjii. "Dirty-boson" systems such as liquid ''He absorbed 
in aerogel, Vycor, or Geltech^, or disordered granular su- 
perconductors^, have also been extensively studied. 

Cold atom gases are becoming important tools to un- 
derstand the interplay between single-particle random- 
ness and cooperative effects such as superfluidity and 
many-body effects induced by interactiongi. Atoms 
trapped in an optical lattice (OL) are particularly suit- 
able candidates for such studies, especially because they 
allow one to reach the strongly-coupling regime through 
the depression of the kinetic energy associated with well- 
to- well tunneling^. A ^^Rb Bose-Einstein condensate 
(BEC) inside a disordered II? OL has been used to study 
the interplay between repulsive interactions and disor- 
derA. In this work it has also been pointed out that 
thermodynamic quantities, such as the superfluid den- 
sity in the case of a BEC, provide a better indicator of 
disorder-induced localization than time-of-flight absorp- 
tion images. The present work has been motivated by 
the experiments in Refs. (Tj andj^ We report a study of 
the interplay between interactions and randomness in a 
repulsive two-component Fermi gas trapped in a II? OL. 
An added motive of interest is that the Landau Fermi- 
liquid paradigm does not apply^" . Two-component Fermi 
gases have recently been prepared in a quasi- II? geome- 
trjiii, thus opening the way to experimental studies of II? 
phenomena such as spin-charge separation-'^ and atomic- 
density wavesi^. 

The ground state of an interacting Fermi gas mov- 
ing under harmonic confinement in a II? OL shows 
in the absence of disorder five qualitatively different 
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How does disorder influence these phases and their ther- 
modynamic properties? In the following we provide a 
quantitative answer to this question. In particular, we 
demonstrate that the incompressible Mott-insulating re- 
gions are very stable against disorder at strong coupling. 
We also show that the compressibility exhibits a disorder- 
induced low-density anomaly, similar in some respects to 
the one which has been found both experimentall}!^ and 
theoretically"^ in the 21? electron liquid close to the onset 
of the MIT. 

The ID random Fermi-Hubbard model — We consider 
a two-component Fermi gas with A^f atoms constrained to 
move under harmonic confinement of strength V2 inside 
a disordered II? OL with unit lattice constant and Ns 
lattice sites i € [IjA^s]- The system is described by a 
single-band Hubbard Hamiltonian, 

H = - Y^^t.j {cl^Cja +}i.c.) + UY^ n,-fhii 
+ 1^2^(«- Ar,/2)2n, + ^e,n,. (1) 

i i 

Here = i > if are nearest sites and zero other- 
wise, a —^tI is a pseudospin-1/2 label for two internal 
hyperfine states, 

^ia — cJ^Cj(j is the pseudospin-resolved 
site occupation operator, and fii = fii^. The effect of 
disorder is simulated by the last term in Eq. J^l, where 
Si is randomly chosen at each site from a uniform distri- 
bution in the range [-1^/2,14^/2] 

In the unconfined limit (V2 = 0) the Hamiltonian TL 
reduces to the Anderson localization problem.^^ for [/ = 
and to the exactly-solvable Lieb-Wu model^fl for = 0. 
The Lieb-Wu model describes a Luttinger liquid away 
from half and full fiUing^^, a Mott insulator at half filling, 
and a band insulator at full filling. In the unconfined 
limit the 2D version of Ti, has been studied in connection 
with the 21? MIT—. In the clean limit the ground-state of 
ii. has been studied by several authorsi^iiSiiSiii, finding 
five different phases as already shown in Fig. 

A particular set ei{a) of randomly chosen values of 
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Si, labeled by a, is a realization of the disorder po- 
tential. Each realization defines an external potential 
Vi{a) = ¥2(1 — iVs/2)^ + £i(a) which determines, via 
site-occupation functional theorjiS, a site occupation 
ni{a) = {'i' alnil'i' a) , I^'q) being the ground state (GS) of 
TC for disorder realization a. The total energy is a unique 
functional of ni{a), which can be written as^iS^ 



£[n{a)] = J^HK[n(a)] 



(2) 



Here .?HK['^(a)] is a universal functional of the site oc- 
cupation, in the sense that it does not depend on the 
external potential Vi{a). The GS site occupation can be 
found by solving the Euler-Lagrange equation 



Sni{a) 



(3) 



the constant fi{a) being a Lagrange multiplier to en- 
force particle conservation, i.e. '^ini{a) = Nf. A 
local-density approximation (LDA) will be used below 
for J^HK [«(«)]• 

The site occupation Mi due to H is finally obtained by 
means of a disorder ensemble average, i.e. Mi = ((ni})dis 

where ((C'))dis = liniM^oo ^Y.^=i^i^)- practice 
one can average only over a finite number M of dis- 
order realizations. Due to our very efficient compu- 
tational method^^ we are easily able to average over 
Ml = 10** realizations of disorder for every set of pa- 
rameters {Ns,Nf,u = U/t,V2/t,W/t}. Averaging over 
such large numbers of realizations becomes necessary as 
the strength of disorder increases and we have checked 
that the density profiles that we report below are sta- 
ble against further increases of A4. Finally, the global 
compressibility can be obtained from the stiffness Sp = 
{{6n/6Nf))dis (see Ref.[l3). 

Site occupation and stiffness anomaly in the presence 
of disorder — The functional J-uk is the sum of three 
terms, 

.^HKKa)] =T,[n(a)]-f |^n2(a)+5,,[n(a)]. (4) 

i 

The first term is the nonintcracting kinetic energy func- 
tional, which is approximated in this work "a la Thomas- 
Fermi"iiSiii, i. e. through an LDA based on the nonintcr- 
acting kinetic energy of the Lieb-Wu model. The other 
terms in J-hk are the mean-field interaction energy and 
the exchange-correlation energy functional incorporating 
many-body effects beyond mean field. This is approx- 
imated through an LDA based on the exactly known 
exchange-correlation energy of the Lieb-Wu modeUl. 
Within this LDA an explicit equation for ni{a) can easily 
be derivediSiii from Eq. Its self-consistent solution 
leads to the site occupation ni{a) and to the chemical 
potential /x(a) for a particular disorder realization^'*, and 
hence to the disorder-ensemble averaged quantities. 



In Fig. 121 we show the site occupation Mi for a disor- 
dered gas with u — 4, which is in phase B ior W — 0. The 
edges of the Mott plateau are the first to be corrupted 
by the appearence of weak disorder. The stability of the 
central region of the plateau depends on the number of 
atoms, if all other parameters are kept fixed. For iVf — 60 
(top panel of Fig.[2Jl the central region persists over a fi- 
nite range of disorder. Disorder increases the tunneling 
through the edges of the trap, making the confinement 
effectively weaker, and thus leads to broadening of the 
site occupation. The Mott plateau has disappeared at 
W/t — 4, and for W/t = 20 the site-occupation profile 
has a strongly non-parabolic overall shape. For iVf = 70 
instead (bottom panel of Fig. 121, the Mott plateau at 
the center of the trap is unstable against the formation 
of a fluid phase with TVi > 1, and the Mott phase can 
survive for weak disorder only in an intermediate region 
between the edge and the center of the trap. In fact, 
N{ = 70 is the critical number of atoms at which the 
phase transition B ^ C occurs in the clean limit: a weak 
disorder potential can shift the transition and induce a 
fluid region embedded in the Mott plateau (see the plot 
corresponding to W/t = 1 in the bottom panel of Fig.[2Jl. 
Eventually, the scenario depicted above for iVf = 60 is 
re-established when W is strong enough. 

In the top panel of Fig. |21 we show the site occupa- 
tion for a disordered gas with w = 8, which is in phase 
V for W ^ Q. Clearly, the band-insulating region can 
only be corrupted from below. We note that at W/t = 3 
the Mott-insulating regions still exist, while the band- 
insulating region has been destroyed. This conflrms the 
expectation that a Mott-insulating region, having its ori- 
gin in exchange-correlation effects, is more stable against 
disorder than a band-insulating region. These behaviors 
are examined in detail in the bottom panel of Fig. |21 

We turn in Fig. 0] to illustrate the effect of disorder on 
the stiffness of the Fermi gas. In the top panel we show 
Sp as a function of TVf at different values of u in the ab- 
sence of disorder. At m > 4 this quantity exhibits three 
non-analyticity points associated with the three phase 
transitions: A ^ B, B ^ C and C — > I?. In the bot- 
tom panel of Fig. 0] we show the same disorder-averaged 
quantity for a disordered Fermi gas with u = 8. We 
see that the disorder has two main effects. It not only 
leads to smoothing of the non-analytic behaviors found 
in the clean limit, but also induces a strong stiffening at 
low density. The latter is an "anomalous" behavior com- 
pared to that found in the clean limit. In fact, for finite 
W the stiffening appears to grow unbounded at very low 
density (see bottom panel of Fig.^J, following the power 
law Sp (X (TVf)"" with an exponent v ~ 0.6. The value 
of the exponent is essentially independent of the param- 
eters u and W/t, but depends on the confinement: for 
example, we find w 0.4 for an open lattice with V2 = 0. 
At high density Sp appears instead to be essentially un- 
affected by disorder. 

The low-density behavior of Sp is reminiscent of what 
has been found in Refs. |3and|3 for a 2D electron liquid. 
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Following Ref. we explain the origin of the anomaly 
using the concept of density percolation. As N{ decreases 
the high-density regions tend to become disconnected, 
since the atoms tend to occupy just the deepest valleys 
in the disorder landscape. At given u and low A^f, the 
system thus stiffens as the disorder grows (see the bottom 
panel in Fig.^. For a given realization of disorder, the 
number JVq of essentially empty sites increases with W/t, 
as it is shown in the inset. 

There is, however, an important conceptual difference 
between the present Fermi gas and the 2D electron liq- 
uid. In the latter the density is also an inverse measure 
of the coupling strength: the stiffness anomaly at low 
density occurs in the strongly correlated regime. In the 
present case the atom number A^f and the interaction 
parameter U are instead independent parameters. The 



anomaly that we observe occurs also in the noninteract- 
ing limit (see the bottom panel of Fig.EJ, demonstrating 
the crucial role of the quenching of percolation in origi- 
nating the anomaly. The interatomic repulsions enhance 
the stiffness at low density in the disordered case just as 
they do in the clean case, in accord with the intuitive 
expectation that a repulsive system is less compressible. 

In summary, we have shown how disorder affects the 
quantum phases of interacting Fermi gases moving under 
harmonic confinement in ID lattices. In particular we 
have seen that Mott insulating regions are quite stable 
against uniformly distributed uncorrelated disorder and 
that the disorder induces an anomalous increase of the 
stiffness at low density from quenching of percolation. 
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FIG. 1: (color online) Sketch of the ground-state site occupation Ui of an interacting Fermi gas in a harmonic trap and a 
clean ID lattice. Phase ^ is a fluid (known as the metallic phase in the case of electrons) with < < 2. In phase B an 
incompressible Mott insulator occupies the bulk of the trap with rii locally locked to 1. In phase C a fluid with 1 < < 2 
is embedded in the Mott plateau. In phase O a band insulator with rii locally locked to 2 is surrounded by fluid edges and 
embedded in the Mott plateau. Finally, in phase £ a band insulator in the bulk of the trap coexists with fluid edges. 




FIG. 2: (color online) Site occupation Mi as a function of site position i for u — A and V2/t — 2.5 x 10 ^ in a lattice with 
A^s = 200 sites. The number of atoms is A'^f = 60 in the top panel and N{ — 70 in the bottom panel. 
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FIG. 3: (color online) Top panel: Site occupation TVi as a function of i for Nf = 200 atoms in a lattice with A'^b = 200 sites, in 
the case w = 8 and V2/t = 2.5 x 10~^. Bottom panel: Number of consecutive sites A/iviott (A/sand) at which \Mi — 1| < 10~^ 
{\Ni — 2| < 10~^), as a function of W/t for the system shown in the top panel. The steps show that each insulating region is 
stable over a finite range of W/t. Muott can be fitted with a linear function over the whole range of W/t (solid line through 
the squares, showing that TViviott = at {W/t)c\ ~ 5). A/Band has a linear behavior up to W/t « 1.8 and beyond can only be 
fitted by a nonlinear function (solid and dashed lines through the dots, showing that A/band = at {W/t)c2 ~ 3 < {W/t)ci). 
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FIG. 4: (color online) Thermodynamic stiffness Sp (in units of t) as a function of iVf for V2/t = 2.5 x 10"^ and A^s = 200 lattice 
sites. Top panel: results for a clean system at various values of u (in the noninteracting case only the phase transition A £ 
can occur). Bottom panel: results for a disordered system at u = 8 and for < W/t < 10. The black triangles report the 
low-density stiffness of a noninteracting system at W/t = 10. The inset shows the number Ao of sites at which ni{a) < 10~* 
in a particular realization of disorder, as a function of A^f for u = 8 and < W/t < 10. 



